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c/3 ■ Abstract. The expansion of laser-irradiated clusters or nanodroplets depends 

strongly on the amount of energy delivered to the electrons and can be controlled 

by using appropriately shaped laser pulses. In this paper, a self-consistent kinetic 

^h . model is used to analyze the transition from quasineutral, hydrodinamic-like expansion 

f> |' regimes to the Coulomb explosion (CE) regime when increasing the ratio between the 

thermal energy of the electrons and the electrostatic energy stored in the cluster. It is 
shown that a suitable double-pump irradiation scheme can produce hybrid expansion 
regimes, wherein a slow hydrodynamic expansion is followed by a fast CE, leading to 
ion overtaking and producing multiple ion flows expanding with different velocities. 
This can be exploited to obtain intracluster fusion reactions in both homonuclear 
deuterium clusters and heteronuclear deuterium-tritium clusters, as also proved by 
^. . three-dimensional molecular-dynamics simulations. 
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1. Introduction 

The interaction of ultraintense lasers with jets of molecular clusters or nanodroplets 
(with typical size in the range 1 — 100 nm and containing 10 2 — 10 8 particles) is a 
central research topic [1], with important applications, such as tabletop nuclear fusion 
for compact neutron sources [2-8], or the laboratory investigations of nucleosynthesis 
reactions, relevant to astrophysical scenarios [9-11]. 

Clustered media can be regarded as sparse distributions of tiny solid targets, a 
peculiar configuration that allows for both a deep penetration of the laser radiation and 
a strong laser-matter coupling with many individual, overdense targets, thus providing 
extremely efficient energy absorption [12]. When hit by an ultraintense laser beam, the 
neutral atoms in a cluster are promptly ionized (cf. Ref. [13] for a detailed analysis 
of the concurring ionization mechanisms in different laser/cluster configurations) and a 
dense "nanoplasma" [1,14] is formed. The free electrons then absorb energy from the 
laser pulse [15] and start expanding, causing the formation of strong electric fields, which 
lead to efficient ion acceleration, as first predicted by Dawson [16]. When the energy 
transferred to the electrons is much smaller than the electrostatic energy stored in the ion 
core, charge separation is localized to regions much smaller than the cluster [14,17-21], 
which then remains quasi- neutral and undergoes a hydrodynamic-like expansion [22-26]; 
in opposite conditions (e.g. with small deuterium clusters exposed to extremely intense 
laser radiation) the electrons suddenly escape from the cluster and the remaining 
bare-ion distribution undergoes a pure Coulomb explosion (CE) [27]. In intermediate 
situations, the expansion dynamics is a mixture of the phenomenology of the two limits, 
with the expansion process being strongly dependent on the self-consistent dynamics of 
ions and trapped electrons [21,28]. When increasing the laser energy, or when lowering 
the cluster size and density, the expansion conditions vary smoothly from quasi- neutral, 
hydrodynamic-like regimes to pure CE regimes, as confirmed by particle-in-cell (PIC) 
simulations [29,30] of the self-consistent laser-cluster interaction and by kinetic or fluid 
modeling of the expansion of finite-size, non-quasi-neutral plasma expansions [21,28,31]. 
Therefore, the expansion regime can be controlled by regulating the amount of energy 
transferred to the electrons [32] , which can be obtained with appropriately shaped laser 
beams. An important example is the irradiation of homonuclear deuterium clusters with 
two sequential laser pulses having different intensities [33]. In this way, one can taylor 
the expansion dynamics so as to induce overrunning between ions and the consequent 
formation of expanding shells containing multiple ion flows. Within such structures 
(here denoted as "shock shells", following the terminology in [27]), the relative velocities 
between deuterium ions from different flows can be sufficiently high for energetic 
collisions and dd fusion reactions to occurr [34]. Since these intracluster reactions occur 
early in the expansion, and before different exploding clusters overlap, they are expected 
to produce a time resolved burst of fusion neutrons before the bulk neutron signal due 
to intercluster reactions [27,34]. In the case of heteronuclear clusters, e.g. deuterium- 
tritium clusters, ion species having different charge-to-mass ratios expand with different 
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velocities making intracluster, interspecies reactions possible also with standard single- 
pulse irradiation [11,35]. 

This paper provides an organic review of the work on controlled expansions of 
clusters and nanoplasmas published in [28,31-34], complemented with novel simulation 
results obtained with a recently developed molecular-dynamics technique [36]. In the 
following, the transition to the CE regime is analyzed with a self-consistent kinetic model 
(section [2]), the concept of shock shell is briefly reviewed (section [3]), and the possibility 
of achieving intracluster nuclear reactions in homonuclear clusters [27, 34] by driving a 
slow expansion followed by a sudden CE [33,34] is explored (section H]). 

2. Kinetic analysis of spherical plasma expansions 

In order to analyze the influence of the electrons on the expansion dynamics, a kinetic 
model for the collisionless expansion of a spherically symmetric nanoplasma has been 
developed [28], based on the assumptions that the electrons are nonrelativistic and 
resorting to the large mass disparity between electrons and ions. The expansion process 
is divided in two stages: an initial charging transient with frozen ions, and the long 
term expansion of both ions and electrons. In the second stage, the expansion dynamics 
is described self-consistently by following the radial motion of the cold ions, while 
accounting for the three-dimensional dynamics of the hot electrons using a sequence 
of "ergodic" equilibrium configurations [31], represented by stationary solutions of 
the Vlasov equation that depend on the total energy e = mv 2 /2 + $(r) only. The 
electron density instantaneously in equilibrium with the electrostatic potential $ can 
be written as n e (r; {$}) = j^pi f p (e) & (r, e; {$}) de where {$} indicates functional 
dependence on $, p is the energy distribution of the electrons, and 8? (r, e; {$}) = 
r 2 (e + e$)2 / j r' [e + e$ (r')] 5 dr' is the probability density of finding an electron 
having total energy e at the radial position r. The self-consistent potential $ satisfies the 
nonlinear Poisson equation V 2 $ = Aire [n e (r; {$}) — n^}. By resorting to the theory of 
adiabatic invariants for time varying Hamiltonians [31,37,38], a closed set of equations 
is then obtained in the form 

ffir- (9<J> 

M w = ~ Ze ^ (1 <" 

r 2 dr \ dr J 

V I CJT' 

m(n) = Mro)-^/ -^ (ic) 

n e = ^Jp (e )^(r,e;{^})de (Id) 

— = -e / — ^(r,e;$)dr (le) 



1 ° f r 2 ^-\ = 4ne (n e - Zn { ) (lb) 



where M and Ze are the ion mass and charge, respectively. Equations ( TTaj) and ( TTej) . 
determining the ion trajectory rj(ro,t) (ro: initial radius) and the evolution of the 
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electron energy e(eo,t) (er> initial energy), are coupled with the nonlinear Poisson 
equation (fT5|) . In Eq. ffldj) . the relation p(e, £)de = p (e )de has been used, where 
Po(eo) is the energy distribution of the electrons resulting from the initial charging 
transient, a fast process that cannot be described by Equations (ITaj) - (fTe|) [31, 32]. 
The self-consistent shape of po(eo) is determined by a virtual charging transient, in 
which the ions stay immobile, while an external potential barrier, initially confining 
the electrons, is gradually moved from R to infinity, with a series of small radial 
displacements (the validity of this procedure has been confirmed by particle-in-cell 
simulations [28] and ad-hoc solutions of the VP model [31]). The initial, nonequilibrium 
electron distribution is assumed to be a Maxwellian, with temperature To, so that the 
whole expansion dynamics is fully determined by the single dimensionless parameter 
T = ZkBTo/ecE = 3A D0 /-Rq ( 6 ce = Ze 2 N Q /R : maximum ion energy attainable from 
the CE of a uniform ion sphere with radius Rq and total charge eiVo; Ado : initial 
electron Debye length), which accounts for both the initial electron temperature and 
the cluster size and density. The positive charge buildup at the ion front, AQ, and 
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Figure 1. (Color online) Charge buildup AQ (left panel) and mean kinetic energy of 
trapped electrons, £, (right panel) as functions of To, after the charging transient with 
immobile ions. Circles refer to the ergodic model, solid lines to the fit laws in the text. 
Dashed lines show the power-law behaviors for To <C 1. 



the mean kinetic energy of the trapped electrons, £ , as functions of T , are displayed 
in Fig. [H Simple fits for AQ and £ are found as AQ = eNo^.eo ( \/6/eT ) and 
£ = §AkT [l - JS.35(1-86T 1/2 )], where F^x) = x/(l + x") 1 ^. For T < 1 (i.e. 
Ado ^ -Ro)j the fits reduce to AQ/(eN ) ~ (6T /e) 1//2 , thus recovering the theoretical 
results for planar expansions [18,39], and to £/(§ &bTo) — 1 — 1.86 T . 

When the initial equilibrium is reached and the ions are allowed to expand, gaining 
kinetic energy, the electrons cool down and the charge buildup decreases (until a ballistic 
regime is reached for both species [40]). Depending on the value of To, such behavior 
strongly affects the ion dynamics and the resulting distribution of ion energy [28,31]: 
for T > 0.5, the spectrum is monotonic as in a CE; for T < 0.5, it exhibits a local 
maximum far from the cutoff energy, thus being qualitatively different from the CE case. 
In this respect, the transition value To = 0.5 can be considered as the lower bound for the 
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Figure 2. (Color online) Asymptotic ion energy spectra for different values of 
To, compared with the theoretical asymptotic spectrum for the uniform CE case, 
3Ro/(2Z 2 e 2 ) (e/ecE) (left panel); cutoff ion energy and location of the maximum 
in the ion energy spectrum as functions of To (right panel). Circles refer to the ergodic 
model, solid lines to the fit laws in the text. The dashed line represents the power-law 
behavior of e m ax for To <C 1. 



validity of the CE approximation. The cutoff energy is fit by e max = T\ .43(2.28 T )ecE 



which, for T < 1, reduces to e r 



2.28 T 6ce; for T < 0.5, the location of the 
maximum is fit by the power-law e pea k = 0.3TQ 9 ecE- These scaling laws, valid for any 
combination of Rq, n , and T , can be useful to interpret experimental data [28,41-43]. 



3. Shock shells in nonuniform Coulomb explosions 



In the absence of electrons, the explosion dynamics of a spherical distribution of cold 
ions is described by Equation (fTqj) . which reduces to Md 2 r l jdt 2 = Z 2 e 2 Ni(r[,t)/rf, 
where N^r, t) is the number of ions enclosed by a sphere of radius r at time t. 

If the initial profile of ion density is uniform and step-like (i.e. n i0 for r ^ Rq), 
the repulsive electric field grows linearly for r < Rq, reaching its maximum at the outer 
boundary, which causes the outer ions to be always faster than the inner ones. Thus, ions 
never overtake each other, N[(r,t) is conserved along the ion trajectories (so is the total 
energy of each ion), and the equation of motion can be integrated analytically [27,44-46], 
yielding, for ions initially at rest, the expansion velocity v- x = [2(r ; — ro)/(3r i )] 1 / 2 roco'pi, 
where u pi = (47rZ 2 e 2 n i0 /M) 1 / 2 is the initial ion plasma frequency. The radial trajectory 
is then given, in implicit form, by [£ (f, - I)} 1 ' 2 + log[^ 1/2 + (£ - l) l > 2 } = (2/3) 1 / 2 cj pi t, 
where £1 = rj/Yn is the expansion factor. In this solution, £; is independent of r , the v — r 
phase-space profile is always a straight line with equation v = [2(£ — l)/(3£ 3 )] 1 / 2 c<jpir, 
and the ion density maintains its step-like form, decreasing in time as U[Q^(t)~ 3 . The 
asymptotic ion energy distribution is 3R /(2Z 2 e 2 ) (c/cce) 1 ^ 2 (cf. Fig. [2]). 

If, more realistically, nonuniform density profiles are considered, the expansion 
features change qualitatively [33,34]: when the initial density is a decreasing function 
of r, the repulsive Coulomb field reaches its maximum within the ion sphere [27,30], 
leading to ion overtaking and to the formation of multiple-flow regions (shock shells) 



Expansion of nanoplasmas and laser-driven nuclear fusion in single exploding clusters^ 

with characteristic multi-branched phase-space profile as in Fig. [3] (since the ion 
trajectories are no longer independent of one another, the above analytical solution 
no longer holds [27,30,44]). The physical interest of these shock shells resides in the 
appearance of multiple flows with large relative velocities within a single exploding 
cluster, which can lead to energetic ion-ion collisions and to intracluster reactions in 
homonuclear clusters [27,33,34]. An effective strategy to produce large-scale shock 
shells in a controlled fashion is combining different expansion regimes so that a slow 
hydrodynamic-like expansion, providing a smoothly decreasing ion-density profile, is 
followed by an abrupt CE. This is achievable with a double-pump irradiation scheme 
as in Fig. [3j wherein a weak pulse (intensity: I\\ duration: ri) is followed by a strong 
pulse (intensity: I 2 > I\\ duration: r 2 ), with suitable time delay At: I\ must be 
sufficiently high to ionize the atoms, creating a nanoplasma, but not so high as to expel 
a significant fraction of the electrons from the cluster, whereas J 2 must be high enough 
to drive a sudden CE. In these conditions, a pronounced shock shell is formed, whose 
features are determined by the key double-pump parameters, namely I\ and At, which 
strongly affects the density profile of the cluster when being hit by the second laser. 
The effectiveness of the technique has been demonstrated [32-34] by resorting to two- 
and three-dimensional PIC simulations performed using the OSIRIS 2.0 framework [47], 
closely matching realistic physical scenarios. 
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Figure 3. (Color online) Sketch of the double-pump irradiation scheme (left panel) 
and example of the typical formation and evolution of a shock shell with three-branched 
profile, the shock shell, in the v — r phase space (right panel): early in the explosion, 
the phase-space profile bends on its right until it becomes multi- valued at time t = t s h 
and an expanding shock shell forms, with inner and outer boundaries r s h and R s h, 
respectively. Energetic collisions can occur between ions belonging to different velocity 
branches. 



4. Intracluster nuclear reactions 



The ability of producing multiple flows with high relative velocities within single 
exploding clusters makes the phenomenon attractive as a possible way to induce 
intracluster fusion reactions in homonuclear deuterium clusters. As pointed out in [27], 
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the rates for intracluster reactions can be significantly higher than those for intercluster 
reactions, because the typical densities within a shock shell (< n i0 ~ 10 22 cm" 3 ) are 
higher than those within the hot plasma filament resulting from the exploded clusters 
(~ 10 19 cm -3 ). However, it is crucial to consider also that the expanding shock-shell 
stays appreciably dense only for a very brief time (~ 10 fs), much shorter than the 
typical disassembly time of the plasma filament (~ 10 — 100 ps): as shown in [34], the 
rates of intracluster reactions exhibit a sharp, time-resolved peak right after the shock- 
shell formation. Once the expansion dynamics is known, the number of intracluster 
fusion reactions is obtained by summing over all possible contributions from collisions 
between ions belonging to different velocity branches in a shock shell (cf. Fig. [3]). The 
number of reactions per unit time and volume, 1Z, is given by 

K = ^ n h {r)n k {r)a (\v h - v k \) \v h -v k \, (2) 

h<k 

where a is the cross section for dd fusion, while nt and v h indicate, respectively, the ion 
density and velocity on the hth branch. The number of intracluster fusion reactions, A/", 
is obtained by integrating TZ over time and space, as 

r+oo rR sh 

N = 4vr / / TZr 2 drdt, (3) 

where r s h, -R s h, and t s h are the shock-shell boundaries and formation time, respectively. 
In [34] , the influence of the double-pump parameters on M was analyzed using a simple 
ID model wherein the laser field of the first pulse gradually strips the initially neutral 
cluster of a part of its electrons according to a cluster barrier suppression ionization 
(CBSI) model (cf. Ref. [13]). This provided preliminary estimates for the optimal 
combinations of delay and intensities that maximize A/", and suggested that, for very 
large clusters (i?o > 100 nm), the intracluster reaction yield can become comparable 
with the intercluster neutron yield, with ~ 10% of the fusion reactions arising from 
intracluster collisions [34]. 

Accurate calculations of the intracluster reaction yields achievable with the double- 
pump irradiation of homonuclear deuterium clusters are currently being performed 
resorting to three-dimensional molecular-dynamics simulations based on the recently 
developed scaled electron and ion dynamics (SEID) technique [36,48], for different 
values of laser intensities, durations, and delays. Preliminary results indicate that, 
with I x in the range 10 15 — 10 17 Wcm~ 2 , T\ in the range 20 — 30 fs, I 2 = 10 20 Wcm~ 2 , 
T2 = 20 fs, and At in the range 110 — 180 fs, the number of intracluster reaction per 
cluster can be as high as M = 4 x 10~ 3 with a cluster size Rq = 110 nm. These 
values are comparable with the estimate in [34]; furthermore, as predicted in [34], it 
is found that the optimal intensities of the first pulse are significantly lower, typically 
by an order of magnitude, than those estimated with the CBSI model. Calculations 
based on the SEID technique are being performed also in the case of heteronuclear 
deuterium-tritium clusters and nanodroplets irradiated by a single laser pulse, for either 
homogenous mixtures or layered targets composed by a deuterium core surrounded by 
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a tritium shell. Owing to the higher value of the deuterium-tritium fusion cross section 
for the energy range considered, higher number of intracluster reactions, on the order 
of M = 0.1 — 1, are obtained with respect to the pure-deuterium case. Detailed results 
of these numerical simulations will be presented in a future publication [48] . 

5. Summary 

The ergodic model has been used to simulate the collisionless expansion of spherical 
nanoplasmas driven by energetic electrons over a wide range of values of plasma size, 
plasma density, and initial electron thermal energy, combined in the single dimensionless 
parameter To: a qualitative change in the shape of the asymptotic energy spectrum of the 
ions when approaching the Coulomb explosion regime has been identified and accurate 
fit laws for the relevant expansion features have been provided. The douple-pump 
irradiation technique to induce the formation of multiple ion flows during the expansion 
of large homonuclear deuterium clusters or nanodroplets has been described, and its 
applicability to produce intracluster dd and dt fusion reactions has been proved resorting 
to three-dimensional numerical simulations performed with a recently developed scaled 
molecular dynamics technique. 
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